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Abstract 

We investigate the evolution of the hard X-ray luminosity of the red galaxy 
population using a large sample of 3316 red galaxies selected over a wide range 
in redshift (0.3 < z < 0.9) from a 1.4 deg 2 region in the Bootes held of the 
NOAO Deep Wide-Field Survey (NDWFS). The red galaxies are early-type, 
bulge-dominated galaxies and are selected to have the same evolution corrected, 
absolute i?-band magnitude distribution as a function of redshift to ensure we 
are tracing the evolution in the X-ray properties of a comparable optical pop¬ 
ulation. Using a stacking analysis of 5-ks Chandra/ACIS observations within 
this held to study the X-ray emission from these red galaxies in three redshift 
bins, we hnd that the mean X-ray luminosity increases as a function of red¬ 
shift. The large mean X-ray luminosity and the hardness of the mean X-ray 
spectrum suggests that the X-ray emission is largely dominated by AGN rather 
than stellar sources. The hardness ratio can be reproduced by either an absorbed 
(Nh ~ 2 x 10 22 cm -2 ) r=1.7 power-law source, consistent with that of a popula¬ 
tion of moderately obscured Seyfert-like AGN, or an unabsorbed T=0.7 source 
suggesting a radiatively inefficient accretion how (e.g., an advection-dominated 
accretion how). We also hnd that the emission from this sample of red galax¬ 
ies constitutes at least 5% of the hard X-ray background. These results suggest 
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a, global decline in the mean AGN activity of normal early-type galaxies from 
z ~ 1 to the present, which indicates that we are witnessing the tailing off of the 
accretion activity onto SMBHs in earlv-type galaxies since the quasar epoch. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: elliptical 
and lenticular, c-D — X-rays — cosmology: observations 


1. Introduction 

The red galaxy population hosts the bulk of the stellar mass at low redsliift and is 
predominantly a bulge-dominated population (e.g., Hogg et al. 2002; Kauffmann et al. 2003). 
It is now known that virtually all local bulge-dominated galaxies host a super-massive black 
hole (SMBH) at their center (e.g., Begelman 2003; Kormendy & Richstone 1995). These 
SMBHs are thought to be relics of an earlier accretion phase onto active galactic nuclei 
(AGN). The recent discovery that the mass of the SMBH is strongly correlated with the 
velocity dispersion and mass of the galactic bulge (Magorrian et al. 1998; Gebhardt et al. 
2000; Ferrarese & Merritt 2000) suggests an intimate relationship between the growth of 
the black hole by accretion processes and the build-up of the galaxy bulge. While a tight 
correlation exists between the two today, the relationship at earlier times is less clear, with 
some observational studies favoring near-lockstep evolution (Heckman et al. 2004) while 
others suggest large variations in the relative bulge/black hole growth rates and mass ratio 
(e.g., Grupe & Mathur 2004). To properly understand this relation, one must have better 
measurements of the accretion history onto the central SMBH. 

Studies with the Hubble Space Telescope show that the hosts of luminous quasars are 
generally massive, early type galaxies (e.g., McLeod & McLeod 2001; Dunlop et al. 2003). 
It is thus likely that many red galaxies may host SMBHs which have formed in an earlier 
epoch of accretion as powerful AGN. It is therefore in this population that one may expect 
low luminosity and/or obscured AGN to reside. 

Hard X-ravs provide us with a relatively unbiased measure of the accretion onto SMBHs 
due to their ability to penetrate all but the most dense columns of obscuring gas. In this 
paper, we use X-rays to investigate the possibility that a significant amount of accretion 
onto SMBHs is still occurring within normal early-type galaxies. We present a preliminary 
analysis of the accretion history of a population of red galaxies selected from the 1.4 deg 2 
initial data release of the NOAO Deep Wide-field Survey (NDWFS; Jannuzi & Dev 1999; 
Jannuzi et al. in prep.; Dev et al. in prep.). Our survey demonstrates the power of stacking 
analyses of X-rays from large, optically selected, samples of galaxies. We assume a ACDM 
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(£1m = 0.3, = 0.7) cosmology and Hubble constant of H 0 = 70 km s x Mpc 1 


2. Data 

2.1. The Red Galaxy Sample 

The optical and near-IR emission from red, early-type, bulge-dominated galaxies is 
dominated by old stellar populations, and the evolution of these systems is well-approximated 
by “passively” evolving models, where the bulk of the stars are formed in a burst at high 
redshift (z > 2) and the population subsequently ages gradually with little additional star 
formation (e.g., Stanford et al. 1998; Kochanek et al. 2000; Rusin et al. 2003). The observed 
optical colors are dominated by a strong 4000 Angstrom break and vary in a relatively well 
understood manner; modelling therefore yields accurate photometric redshifts (e.g., Brown 
et al. 2003; Csabai et al. 2003). The expected low levels of star formation imply a low 
contribution to the total X-ray luminosity from sources such as high-mass X-ray binaries 
which trace the star formation rate (SFR). 

We selected our sample from a 1.4 deg 2 area within the NDWFS Bootes field. The 
NDWFS comprises deep optical and infra-red imaging in two 9 deg 2 regions of the sky, one 
in Bootes and one in Cetus, and the area considered in this paper lies in the northwest 
corner of the Bootes field. The NDWFS is an ideal survey for identifying large homogeneous 
populations due to its deep coverage over large areas. 

Our red galaxy sample selection is described in detail by Brown et al. (2003). We used 
Source Extractor, version 2.2.2 (Bertin & Arnouts 1996) to detect the galaxies and measure 
their optical magnitudes. We then estimated photometric redshifts for the sample using a 
customized routine which simulates the BwRI color-color space using PEGASE2 spectral 
evolutionary synthesis models (Fioc & Rocca-Volmerange 1997) that assume exponentially 
declining star formation rates with e-folding times between r = 0.6 and 15 Gvr (“r- models”), 
and a formation redshift of Zf ~ 4. Comparison with spectroscopic redshifts demonstrate 
that the photometric redshift estimates for this sample are reliable out to z ~ 1 (they have 
a 5% systematic error and a ±5% la uncertainty; Brown et al. 2003). For the red galaxy 
sample, we chose galaxies with SEDs well-fit by models with r < 4.5 Gyr, to select objects 
which map the red envelope of the galaxy color-color distribution and are likely to trace the 
evolutionary history of modern-day ellipticals. The rest-frame (z=0) color of the r = 4.5 Gyr 
model with E(B — V) 0.04 intrinsic dust extinction is B\y — R = 1.44. This is consistent 
with the Sab template of Fukugita et al. (1995); therefore our r < 4.5 Gyr criteria roughly 
selects galaxy types Sab and earlier. 
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We estimated the evolution corrected present-day f?-band absolute magnitude, Mr, of 
each of the galaxies. The PEGASE2 r = 1 Gyr models and observations (Jorgensen et al. 
1999; Schade et al. 1999; Im et al. 2002) suggest the magnitude correction is relatively small: 
z=0.9 red galaxies are ~ 0.9 magnitudes brighter in rest-frame J?-band than z =0 red galaxies. 
We then imposed an absolute magnitude cut Mr < —21.3 to ensure that the faintest galaxies 
in the lowest redshift bin were detectable out to a redshift of z= 0.9, the limit of our sample. 


2.2. X-ray Imaging: The XBootes Survey 

The X-rav data used in this paper are from the Chandra observations of the Bootes 
field of the NDWFS (Murray et ah, submitted; Renter et al. submitted). A 9 cleg 2 region of 
sky, designed to match the area covered with NDWFS Bootes field, was observed in the 0.5- 
7.0 keV energy range by the Advanced CCD Imaging Spectrometer (ACIS) on the Chandra 
X-ray Observatory. This comprised 127 pointings of 5 kilosec each. The X-ray number 
counts and angular correlation function of the detected sources are reported in Renter et 
al. (submitted). Brand et al. (in prep.) presents the optical counterparts to the X-ray 
detected sources. Here we analyze only 1.4 deg 2 (roughly 21 ACIS pointings) of the survey 
area corresponding to the northwest corner of the NDWFS Bootes field, roughly the area 
covered by the DR2 release of the NDWFS (http://www.noao.edu/noao/noaodeep). 


3. Stacking Analysis 

We investigated the evolution in the mean properties of the red galaxies as a function 
of redshift using a stacking technique (e.g., Brandt et al. 2001b; Nandra et al. 2002). This 
technique is based on the assumption that a large sample can provide an ensemble average 
measurement of the flux from a typical member even if any single galaxy in a given population 
is undetected by a particular observation. In other words, one can increase the effective 
observation depth per object by a factor of the sample size: a sample of 10 3 galaxies effectively 
increases the depth of our 5 kilosec Chandra exposure to that of a 5 Msec observation on the 
mean object. In addition, by stacking the X-ray emission from a large sample of optically 
selected host galaxies, one obtains an averaged signal which will not be affected by AGN 
variability of individual objects. 

We divided the sample of 3316 red galaxies into three redshift bins (0.3 < z < 0.5, 
0.5 < z < 0.7, and 0.7 < z < 0.9, consisting of 1040, 1271, and 1005 galaxies respectively) 
using the photometric redshift estimates described above. Figure 1 shows the sum of the 
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redshift likelihood distributions for the galaxies in each redshift range. Although the errors 
in the photometric redshifts results in some cross-talk between adjacent redshift bins, they 
are accurate enough that the lowest and highest redshift bins should be truly independent. 
The mean present-day rest-frame R-band absolute magnitude is constant (Mr=— 20.4) for all 
redshift bins. We created a stacked X-ray image for each redshift bin by centering on each of 
the optical galaxies in the subsample, effectively extracting the appropriate sub-image from 
the Chandra image, and summing all the sub-images. 

When considering a stacked signal, it is important to ensure that the measurements are 
not dominated by individual objects. To check this, we identified the red galaxies within 5 
arcsec of an X-ray detection above the total-band (0.5-7.0 keV) false-probability threshold 
of 5x10 5 and >4 normalized counts (Murray et al. submitted; Renter et al. submitted). 
Assuming a standard power-law model spectrum with photon index T=1.7, this detection 
limit corresponds to luminosities of 3.4 xlO 42 ergs s _1 , 8.6 xlO 42 ergs s _1 , and 16.7 xlO 42 
ergs s” 1 in the 0.3 < z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9 redshift bins respectively. 
There are 8, 15, and 11 X-ray detected red galaxies, which contribute 45%, 49%, and 44% 
of the total signal in the 0.3 < z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9 redshift bins. 
Excluding these sources decreases the signal in approximately the same proportion in all 
redshift bins, roughly preserving the hardness ratio 1 (see Table 1). Although our analysis 
described below excluded the detected sources, our conclusions are only strengthened by 
their inclusion. 

After removing all detected sources, Figure 2 shows that there is a significant detection 
in all three redshift bins (S/N = 9.7, 4.1, and 5.5 in each of the 0.3 < z < 0.5, 0.5 < z < 0.7, 
and 0.7 < z < 0.9 redshift bins assuming a 1.5 arcsec aperture radius). Radial profiles of 
the data are shown in Figure 3. Our measurements are presented in Table 1. 

We measured the X-ray counts from each stacked image using a centered circular aper¬ 
ture with a fixed angular radius of 1.5 arcsec. We estimated the background by performing 
aperture photometry on the stacked X-ray image at a angular radius between 10 and 30 
arcsec from each red galaxy. After excluding all galaxies that were individually detected in 
the X-ray, we subtracted our estimate of the background flux from that of the remaining 
stacked sample. We converted the net flux to an X-rav luminosity using the central redshift 
of the redshift bin. 

Ideally, one would use an aperture radius which has the same physical radius on the 
sky at different redshifts. This would ensure that the X-ray signal was not increasingly 


1 HR = (Ch — C s )/(Cft+C s ) where C h and C s are the counts in the hard (2-7 keV) and soft (0.5-2 keV) 
bands respectively. 



6 



Fig. 1.— The sum of the redshift likelihood distributions for all red galaxies within the 0.3 < 
z < 0.5 (solid line), 0.5 < z < 0.7 (dashed line), and 0.7 < z < 0.9 (dotted line) redshift bins. 
The redshift distributions of adjoining photometric redshift bins have significant overlaps, 
but the uncertainty is small enough that the lowest and highest redshift bins should be truly 
independent. 
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0.50 < 2 < 0.70 


0.70 < 2 < 0.90 



Fig. 2.— The stacked full-band (0.5-7 keV) X-ray images of the individually undetected red 
galaxies with photometric redshifts 0.3< z <0.5 (left), 0.5< z <0.7 (center), and 0.7< z <0.9 
(right). The images are smoothed by a box-car average of 3 pixels. Overplotted is the 1.5 
arcsec aperture radius used to extract the signal (see text). 
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Fig. 3.— Radial surface brightness plots (counts/arcsec 2 ; left) and plots showing the frac¬ 
tional encircled energy as a function of radius (right) for both the data (stars) and integrated 
PSF models (heavy solid line). These plots are shown for the stacked red galaxies with pho¬ 
tometric redshifts 0.3< z <0.5 (top), 0.5< z <0.7 (middle), and 0.7< z <0.9 (bottom). 
The data was background subtracted using aperture photometry. For the surface brightness 
profile, the PSF model was normalized to have the same number of total counts as found in 
the data out to a radius of 5 arcsec. For the fractional encircled energy plots, the data was 
normalized to have the same fractional encircled energy as the model at 5 arcsec. All errors 
shown are Poisson counting errors. The vertical dotted line shows the 1.5 arcsec aperture 
radius in which we performed the aperture photometry. 
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Table 1. Table showing the total X-ray counts, the estimated background X-ray counts, 
the resultant X-ray counts, the detection significance (a), and the average X-ray 
Luminosity in the total (t; 0.5-7 keV), soft (s; 0.5-2 keV), and hard (h; 2-7 keV) bands for 
each redshift bin. The X-ray luminosity was calculated assuming a power-law spectrum 
with photon index T=1.7 (PL Lum), a power-law spectrum with photon index T=1.7 and 
an HI column density of 2xl0 22 cm -2 at the mean redshift of the bin (Abs Lum), and a 
power-law spectrum with photon index T=0.7 (r=0.7 Lum). All models include absorption 
the Galactic HI column density of 1.75xl0 20 cnr 2 . Also shown is the mean total (0.3-8 
keV) band X-ray luminosity predicted for LMXBs (Kim et al. 2004) and the mean hard 
(2-8 keV) band X-ray luminosity predicted for HMXBs (Grimm et al. 2003). 


Redshift 

0.3-0.5 

0.5-0.7 

0.7-0.9 


All sources 

Number of galaxies 
mean R magnitude 
mean Mr 

Counts (t/s/h) 

Bkgd counts (t/s/h) 

Counts-bkgd (t/s/h) 

Detection Significance (cr) (t/s/h) 
PL Lum (10 41 erg/s) 

Abs Lum (10 41 erg/s) 
r=0.7 Lum (10 41 erg/s) 

1040 

19.75 

-21.2 

1424:11.9/734=8.5/694=8.3 

26.74=5.2/10.44=1.8/15.94=4.0 

115.34=13.0/62.64=9.1/53.14:9.2 

22.3/19.4/13.2 

2.34=0.2/0.74=0.09/2.44=0.4 

3.44=0.3/0.34=0.04/1.34=0.2 

2.64=0.3/0.54=0.06/2.04=0.3 

1271 

20.95 

-21.2 

1104:10.5/474:6.9/634=7.9 

33.14=5.8/12.74=3.6/20.34=4.5 

76.94=10.5/34.44=7.7/42.74=9.1 

13.4/9.7/9.5 

3.54=0.5/0.94=0.2/4.54=0.8 

4.84=0.7/0.44=0.1/2.14=0.4 

3,5±0.5/0.54=0.1/3.34=0.6 

1005 

21.85 

-21.2 

974=9.8/524=7.2/454=6.7 

26.64=5.2/10.34:3.2/16.0±4.0 

70.44=11.1/41.84=7.9/29.04=7.8 

13.6/13.0/7.3 

7.34=1.0/2.44=0.4/7.04=1.6 

9.44=1.3/1.34=0.2/3.04=0.7 

6.64=0.9/1.34=0.2/4.54=1.1 

X-ray non-detected sources 

Number of galaxies 
mean R magnitude 
mean Mr 

Counts (t/s/h) 

Bkgd counts (t/s/h) 

Counts-bkgd (t/s/h) 

Detection Significance (<r)(t/s/h) 
PL Lum (10 41 erg/s) 

Abs Lum (10 41 erg/s) 

T=0.7 Lum (10 41 erg/s) 

1032 

19.75 

-21.2 

76 r 8.7 37.r:fi. 1 /39 : 0.2 

26. 1:1:5.1 10. 1 = 3.2 15.8 ! 1.0 
49.64=10.1/26.7d=6.9/23.24=7.4 

9.7/8.3/5.8 

1.04=0.2/0.34=0.1/1.14:0.3 

1.54:0.3/0.14:0.03/0.64=0.2 

1.14=0.2/0.24=0.05/0.94=0.2 

1256 

20.95 

-21.2 

56 ! 7.5 26:5.1 304:5.5 

32.74=5.7/12.54=3.5/20.0±4,5 
23.34=9.4/13.54=6.2/10.04=7.1 

4.1/3.8/2.2 

1.14=0.3/0.34=0.1/1.14=0.6 
1.54=0.5/0.24=0.07/0.54=0.3 
1.14=0.3/0.24=0.08/0.8±0.4 

994 

21.85 

-21.2 

544=7.3/234=4.8/314=5.6 

26.14=5.1/10.14=3.2/15.84=4.0 

27.94=8.9/12.94=5.8/15.24=6.8 

5.5/4.1/3.8 

2.94=0.8/0.84=0.3/3.74=1.4 

3.84=1.0/0.44=0.2/1.64=0.6 

2.74=0.7/0.44=0.15/2.44=0.9 

LMXB Lum (10 41 erg/s) 

HMXB Lum (10 41 erg/s) 

0.174=0.07 

0.10 

0.174=0.07 

0.14 

0.184=0.07 

0.26 
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dominated by off-nuclear radiation such as diffuse emission from hot gas and X-ray binaries 
in more distant sources. However, we chose an aperture of fixed angular size because the 
effects of losing signal outside the aperture radius due to the Chandra point spread function 
(PSF) dominate over a contribution from extended emission and a variable size aperture 
would degrade the signal-to-noise (S/N) at higher redshifts. To justify the choice of a fixed 
angular aperture radius of 1.5 arcsec, we ran Monte Carlo simulations to ensure any input 
redshift evolution in luminosity is correctly recovered using this aperture. These tests are 
described in more detail in Section 4.2. Performing the analysis with a fixed physical radius 
does not significantly change our results. 


3.1. Correcting for the Chandra PSF 

The Chandra point spread function (PSF) broadens markedly as a function of off-axis 
angle (the 50% encircled energy radius increases from 0.5 arcsec on-axis to 6.5 arcsec at an 
off-axis angle of 10 arcmin). This will result in an increasing fraction of the total X-ray 
counts falling outside the aperture radius for galaxies which happen to lie at a larger off-axis 
angle from the Chandra pointing center. We corrected for this effect in the following way. 
Based on a fit to the numerical rav-trace calculations presented in the Chandra proposers’ 
observatory guide, we adopted a model of the Chandra 50% encircled energy radius, r 5 o at 
a given off-axis distance, A 6: 


r 50 = 0.43 - O.1A0 + O.O5A0 2 (1) 

where r 50 is in arcsec and A 9 is in arcmin. We then calculated the PSF for each of the red 
galaxies (each with a different off-axis distance), summed the individual PSFs to create the 
effective PSF for each stacked image, and then integrated to find the fractional encircled 
energy as a function of radius. Figure 3 shows that the shape of the PSF model we used 
provides a good fit to the data, suggesting that our assumption that the X-ray sources are 
point sources is a reasonable approximation. The fractional 50% encircled energy within a 
1.5 arcsec radius, / 50 are 0.415, 0.371, and 0.398 for the 0.3 < z < 0.5, 0.5 < z < 0.7, and 
0.7 < z < 0.9 redshift bins respectively (i.e., we detect only approximately 40% of the total 
flux due to the effects of the PSF). To account for the losses due to the PSF, we applied the 
correction factors of I// 50 , to our flux measurements in each bin. 
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3.2. X-ray Count Distribution 

Figure 4 shows the distribution of counts contributing to the stacked images and iden¬ 
tifies the detected sources. After subtracting the individually detected sources, only ~ 7%, 
4%, and 5% of the galaxies contribute to the counts in our stacked signal in the 0.3 < z < 0.5, 
0.5 < z < 0.7, and 0.7 < z < 0.9 reclshift bins respectively. However, assuming all sources 
have the average PSF, there may be as many as 17%, 11%, and 13% which would contribute 
to the signal if all their counts didn’t fall outside the 1.5 arcsec aperture radius. What we 
cannot tell from the counts is the nature of the true luminosity distribution. The stacking 
analysis by construction assumes that the one- and two-count sources are merely the bright 
tail in the Poisson distribution of the X-ray flux, that the detected photons are drawn at 
random from the red galaxy population, and that the best estimate of the mean luminosity 
is derived by dividing the detected flux among all the stacked sources. This assumption 
is supported by the optical properties of the individually detected and undetected sources. 
We found no significant differences when comparing the distribution of optical absolute 
magnitudes and colors among the sources contributing one or more counts and the overall 
population. As far as we can tell, given the current sample, the individually detected and 
undetected sources appear to be drawn from the same optical population. If only a fraction 
of the red galaxies contributed to the stacked signal, then the mean flux of the X-rav-bright 
sub-population would be larger. Distributing the flux equally over the entire population, as 
we have done, therefore places a conservative lower bound on the mean flux. 

To investigate the possible distribution of sources below our detection limit, some insight 
can be gained by considering the individually detected sources in the Chandra Deep Fields 
(Brandt et al. 2001a; Giacconi et al. 2002). Szokoly et ah (2004) present the spectroscopic 
follow-up of the X-ray sources detected in the 942-ks exposure of the Chandra Deep Field 
South (CDF-S). The CDF-S is a factor of ~ 5 shallower than our stacked signal and ~ 25 
times smaller in area, but has the advantage that the properties of the individually detected 
sources can be distinguished. The CDF-S sample contains 15 X-ray-detected galaxies in the 
reclshift range 0.3 < z < 0.9 whose optical spectra show no obvious signs of AGN activity 
and resemble those of our red galaxy population. The combined hardness ratio of the CDF-S 
sources is 0.07, which suggests obscured AGN activity (e.g., Norman et al. 2004). This is 
consistent with the mean hardness ratios of —0.03, 0.15, and —0.07 for the red galaxies in 
the 0.3 < z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9 reclshift bins respectively. Assuming a 
power-law spectrum with photon index r=1.7, the mean luminosity in the 0.5-7 keV band 
of the CDF-S sources is ~ 1 x 10 42 ergs s _1 . This corresponds to 0.3 counts in our 5-ks 
exposure. 

All 15 detected CDF-S sources (Szokoly et al. 2004) would fall below the 4-count X-ray 
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Fig. 4.— Histograms showing the number of photons detected within the 1.5 arcsec aperture 
radius as a function of the number of red galaxies for the 0.3< z <0.5 (left), 0.5< z <0.7 
(center), and 0.7< z <0.9 (right) redshift bins. Filled bins show the distribution of the 
sources which are individually detected in the X-ray (>4 total counts using WAVDETECT). 
These sources may have less than 4 counts within the 1.5 arcsec aperture radius due to PSF 
effects. The 10 count bin includes all sources with >10 counts (there is one detected source 
in the 0.3< z <0.5 redshift bin with 24 X-ray counts). The dotted lines show the expected 
number of background counts. 
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detection limit of our XBootes survey data. However, the count distribution of these 15 
sources can be used to predict the number of 1, 2, and 3 photon count sources we would 
expect in our 5 ks observation. Assuming that our sample has exactly the same X-ray flux 
distribution as the 15 CDF-S sources, and correcting for our 1.5 arcsec aperture radius, 
we found that the number of sources anticipated from the CDF-S sample is similar to the 
number observed. The CDF-S data are therefore consistent with the hypothesis that the 
X-ray flux in our sample arises from a large fraction of the red galaxies. However, these 
numbers are presently very uncertain because of the small number of photons in our stacked 
data and the limited public studies available from the deep surveys. 

Analysis of the full (9 deg 2 ) survey will greatly increase our ability to study the pre¬ 
dicted photon count distribution. When combined with the growing body of spectroscopic 
information from the Chandra deep surveys, these may place significant constraints on the 
underlying source distribution. We therefore defer a more detailed analysis to a future paper. 
Here, we simply state that the constraints we can currently place from the published deeper 
Chandra surveys are consistent with our assumption that the stacked signal is representative 
of the overall red galaxy population, and not simply the result of a small, X-ray luminous, 
outlier population. 


4. Redshift Evolution of the Mean X-ray Luminosity 

We can use our three sub-samples to investigate the evolution in the mean X-ray lumi¬ 
nosity of the red galaxy population. Once we have determined this, we would like to test 
that no systematic effects are introducted by our stacking method. 


4.1. The Observed Redshift Evolution 

In order to convert the X-ray counts to an X-ray flux, one must assume a model for 
the spectral energy distribution. When we assume a canonical AGN power-law spectrum 
with photon index r=1.7 and the Galactic HI column density of 1.75xl0 20 cm -2 (Stark 
et al. 1992), we significantly underestimate the fraction of photons in the hard-band relative 
to the total-band (0.25 versus the observed fraction of 0.5). This is clearly inconsistent 
with our data and indicates that the shape of the assumed spectral energy distribution is 
incorrect. If we instead assume an intrinsic power-law spectrum with photon index r=1.7, 
the relative count rates in the hard and soft bins can be explained in a consistent scenario 
whereby the continuum is absorbed by an HI column density of ~ 2 x 10 22 cm” 2 (which will 



13 


absorb a much larger number of the soft photons) at the redshift of the galaxy in all redshift 
bins. Alternatively, the relative count rates can be explained by an unabsorbed power-law 
spectrum with photon index T=0.7. 

Figure 5 shows the evolution of the mean X-ray luminosity of the red galaxy population 
as a function of redshift in the total, soft, and hard X-ray bands, based on these two spectral 
models. Both the total and hard band luminosities increase significantly with redshift in 
both cases. Our conclusions are not strongly dependent on the assumed spectral model. If 
we include the contribution from the detected sources, the overall trend does not change (cf., 
Table 1). 


4.2. Monte Carlo Simulations 

In order to assess our ability to recover the true redshift evolution of the X-ray luminosity 
evolution using the stacking approach, we performed Monte Carlo simulations of our stacking- 
analysis. We generated synthetic red galaxy catalogs distributed randomly within an ACIS 
field. We assumed a simple input model for the evolution of total X-ray luminosity with 
redshift (L x oc (1 + ^) 3 ), normalized to the observed mean luminosity in the 0.3 < z < 0.5 
redshift bin. This is a trend that we might expect from Barger et al. (2001). The X-ray 
photons from each galaxy were distributed using a Gaussian PSF model that varies in its 50% 
encircled energy radius as a function of off-axis distance according to Eqn. 1. We modeled 
the ACIS background by distributing photons at random positions in the field based on the 
observed background surface densities of 0.0032 total (0.5-7 keV) photons per square arcsec. 

Having simulated the galaxy and photon catalogs, we then followed the same stacking- 
procedure that was employed for the real data. Figure 6 shows the assumed and recovered 
luminosity evolution based on a 1.5 arcsec radius aperture (r ap ) and demonstrates that we 
successfully recover the input luminosity profile once the PSF correction factors (discussed 
in Section 3.1) have been applied. Although, adopting larger values of r ap enabled us to 
recover a larger fraction of the signal, this introduced greater background contamination and 
correspondingly larger measurement errors. Because we can estimate the PSF corrections 
with reasonable accuracy, we conclude that the combined use of a smaller aperture and PSF 
corrections provides the best strategy for stacking the emission of point-source dominated 
galaxies, like those in our sample. 
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Fig. 5.— Redshift evolution of the mean X-ray luminosity (excluding individually X-ray 
detected sources) and the associated Poisson counting uncertainties for two different spectral 
energy distributions. The X-ray luminosity was calculated assuming a power-law spectrum 
with photon index r=1.7 and an additional HI column density of 2xl0 22 cm -2 at the mean 
redshift of the bin (left) and with an alternative photon index, r=0.7 (right). In both 
cases, we assumed a Galactic HI column density of 1.75xlO 20 cm -2 . The total, hard, and 
soft luminosities are represented by the solid, dashed, and dotted lines respectively. The 
luminosities obtained when including the detected sources can be found in Table 1. 
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Fig. 6.— The redshift evolution of the mean X-ray luminosity assumed in the Monte Carlo 
simulations (solid line) and recovered using using a 1.5 arcsec aperture radius both with 
(dashed line) and without (dotted line) applying the PSF corrections discussed in Section 3.1. 
The error bars in the plot were calculated using bootstrap resampling. 

5. Discussion 

5.1. X-ray Emission Mechanisms 

The total X-ray emission from galaxies may contain contributions from a variety of 
sources: stellar objects, such as low-mass and high-mass X-ray binaries (hereinafter LMXBs 
and HMXBs); diffuse hot gas; and AGN (e.g., Fabbiano 1989). The LMXB population is 
long-lived and generally traces the old stellar population in galaxies and hence the stellar 
mass of the galaxy; in contrast, the HMXBs have short lifetimes and the luminosity of this 
population reflects the star-formation rate in galaxies. 


5.1.1. Contribution of Stellar Sources 

We can estimate the integrated X-ray luminosity expected from LMXBs from relations 
derived from high-spatial-resolution X-ray observations of nearby galaxies. Adopting the 
relation given by Kim et al. (2004) between the absolute K magnitude of early-type galaxies 
and their integrated LMXB X-ray luminosity, and assuming that this relation does not vary 
with redshift, we predict that LMXBs contribute a mean X-ray luminosity of L LMX b ~0.2 
xlO 41 ergs s' 1 in all redshift bins (see Table 1). We note that using the relation of Col¬ 
bert et al. (2004) results in luminosities a factor of ~10 smaller but this is probably an 
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underestimate because they only calculate the integrated luminosity of discrete sources. A 
comparison with Table 1 and Figure 5 illustrates that this contribution is small compared 
with the mean X-ray luminosity from the red galaxy sample; summing over all the galaxies, 
we find that LMXBs should contribute ~ 11%, 9%, and 4% of the stacked total luminosity 
in the 0.3 < z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9 reclshift bins respectively. 

If the number of LMXBs within a galaxy is larger at z = 0.3 — 0.9 than in the local 
Universe, their contributions to the total stacked luminosity could be higher. White & 
Ghosh (1998) suggest that the X-ray luminosity of galaxies at z ~0.5-l should be an order 
of magnitude higher than in the local Universe due to the effects of the declining cosmic- 
star formation rate on the evolution of the LMXB population. However, stacking analyses 
of optically bright HDF-N galaxies at z ~0.5 show that the X-ray luminosities of these 
galaxies can only be up to a factor of ~2 times higher than local galaxies with similar 
optical properties (Brandt et ah 2001c) and models have been revised accordingly (Ghosh & 
White 2001; Ptak et al. 2001). Including such an evolution suggests that the contribution 
of LMXBs to the total stacked luminosity is at most ~ 20%. 

We can also estimate the contribution to the X-ray luminosity from HMXBs by scaling 
from the star-formation rate (SFR) estimates derived from the population synthesis model 
fits to the optical photometry. The mean SFRs estimated for the 0.3 < z < 0.5, 0.5 < z < 0.7, 
and 0.7 < z < 0.9 redshift bins are 1.8 M 0 yr _1 , 2.3 M 0 yr _1 , and 4.1 M 0 yr -1 . For each 
galaxy, we used the scaling of Grimm et al. (2003) in either the high-SFR (> 4.5 M 0 yr _1 ; 
linear) or low-SFR (< 4.5 M 0 yr _1 ; non-linear) regime to estimate the hard X-ray luminosity 
from HMXBs. We calculate hard X-ray luminosities of L H mxb ~ 0.10 xlO 41 ergs s” 1 , 0.14 
xlO 41 ergs s -1 , and 0.26 xlO 41 ergs s _1 originating from HMXBs: i.e., ~ 7%, 9%, and 7% 
of the stacked hard X-ray luminosity in each of the three redshift bins. 

There may also be a small contribution in the X-ray emission arising from hot gas that 
typically surrounds elliptical galaxies (Forman et al. 1994). This emission is predominantly 
in the soft band and thus subtracting any contamination would make our underlying X-ray 
spectrum even harder. We note, however, that a comparison of the radial profiles with the 
PSF shown in Figure 3 show no compelling evidence for a significant extended component. 

We conclude that the large mean X-ray luminosity (~10 41 ergs s 1 ) and hard spectrum of 
the red galaxies strongly suggests that the signal must be dominated by X-rav emission from 
accretion onto a central AGN rather than X-ray emission originating from stellar sources. 



5.1.2. AGN Accretion Processes 


We have shown in Section 4.1 that the mean red galaxy hardness ratio can be reproduced 
by either an absorbed (Nh ~ 2 x 10 22 cnr 2 ) T=1.7 power-law source or an unabsorbed source 
with a harder (T=0.7) spectrum. 

An absorbed (N H ~ 2 x 10 22 cm -2 ) power-law spectrum with a canonical slope, T=1.7 
is consistent with that of a population of moderately obscured, Seyfert-like AGN (Nandra 
& Pounds 1994; Risaliti et al. 1999), whose primary emission mechanism is thought to be 
inverse-Compton scattering originating from a hot accretion disk. The optical and soft X- 
rav photons are thought to be absorbed by an intervening column of gas and dust resulting 
in a hard X-ray spectrum and no optical AGN signature. Ueda et al. (2003) find that X- 
ray absorbed AGNs make up a larger fraction of fainter AGN populations, suggesting that 
absorption may indeed be significant in our sample. 

A harder (T=0.7) spectrum could be produced by accretion onto a SMBH via a ra- 
diatively inefficient accretion flow, such as an advection-dominated accretion flow (ADAF), 
which has a harder characteristic X-ray spectrum than other mechanisms and thus requires 
little or no obscuration to account for the hard spectrum. ADAFs were originally proposed 
by Ichimaru (1977) and Rees et al. (1982), and have been invoked to explain the lack of 
strong X-ray emission expected from accretion of the extensive hot gas known to exist in 
elliptical galaxies (e.g., Naravan & Yi 1995; Quataert & Narayan 1999; Blandford & Begel- 
man 2004). Allen et al. (2000) and di Matteo et al. (2000) analyzed the X-ray spectra of 
6 nearby giant elliptical galaxies observed with ASCA and found a best-fit photon index of 
r=0.6-1.5 and suggest that this is best explained by a ADAF. However, higher resolution 
imaging by Chandra has shown that much of the hard X-ray emission spatially unresolved 
by ASCA originates from luminous LMXBs (Angelini et al. 2001; Loewenstein et al. 2001). 

Although the true underlying X-ray luminosity distribution of the red galaxy population 
is unclear, it appears likely that the majority of the population emits a hard X-ray spectrum 
at a luminosity that is too high to be due to stellar sources and is instead powered by accretion 
onto a SMBH. The results of several studies suggest that only half of all local galaxies have 
evidence for nuclear point sources and/or exhibit AGN signatures in their optical spectra 
(e.g., Ho et al. 1997; Sarajedini et al. 1999; Ravindranath et al. 2001). However, this fraction 
is uncertain due to difficulties in resolving nuclear point sources, distinguishing them from 
compact nuclear star clusters, and identifying AGN in spectra dominated by the host galaxy 
light. In addition, if the accretion flow has become advection dominated, the radiative 
signature is likely to be too weak to detect at optical wavelengths. 
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5.2. Contribution to the Hard X-ray Background 

If obscured AGN and/or low radiative efficiency accretion flows onto SMBHs exist in 
most early type galaxies, their integrated emission may contribute significantly to the hard 
X-rav background (di Matteo & Allen 1999). 

Approximately 94% and 89% of the soft and hard cosmic X-ray background (CXB) 
respectively has now been resolved into discrete sources (Moretti et al. 2003). The CXB 
spectrum is well described by a photon index of T ~1.4 (e.g., Mushotzky et al. 2000). To 
fit the shape of the CXB, models require an as yet unresolved population of hard spectrum 
sources; typically, a population of heavily obscured AGN is invoked (e.g., Ueda et al. 2003). 

The total hard X-ray flux of our red galaxy population (including X-ray detected sources) 
is ~ 1 xlO -12 ergs s -1 cm -2 deg -2 . This sample therefore contributes approximately 5% 
of the hard CXB. This is a lower limit as it only includes a sub-sample of normal galaxies 
optically dominated by stellar light within the redshift range 0.3< z <0.9. More luminous 
AGN within similar galaxies will have been excluded from our analysis because they signif¬ 
icantly modify the optical colors of their hosts and cause them to drop out of our sample. 
The fraction of the CXB contributed by all bulge dominated galaxies is therefore expected 
to be much higher. 

The majority of our red galaxy sample would be detected at the depth of the Chandra 
Deep fields. However, if a large fraction of normal galaxies harbor very low luminosity 
AGN (i.e., below the limit of the deep surveys), they could contribute non-negligiblv to the 
unresolved portion of the hard CXB, because of their high number density. 


5.3. The Nuclear Accretion History of Red Galaxies 

Our stacking analysis of the red galaxy population has demonstrated that their mean 
X-ray luminosity decreases with time. This is perhaps as expected: the dramatic fall in 
the space density of powerful AGN is well known, and has been attributed to a decrease in 
the fuel supply, fueling rate, and/or radiative efficiency onto the SMBH (e.g., Cavaliere & 
Vittorini 2000). We observe a redshift evolution in the total X-ray luminosity of (l + ^) 4 ' 0±2 4 . 
This is consistent with the ~ (1 + z ) 3 variation in accretion activity inferred by Barger et al. 
(2001) from X-ray sources in the CDF-N, and, in corroboration with the above observations, 
is indicative of a significant decline in the mean rate of accretion onto SMBHs in early-type 
galaxies from z ^ 1 to the present. 

However, we need to be careful when comparing our results to those derived from X-ray 
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selected samples. Recent results show that the peak in the number density of luminous 
AGN occurs at higher redshift than that of the less powerful AGN (Cowie et al. 2003; Ueda 
et al. 2003): the so-called “cosmic down-sizing” effect. This suggests that at lower redshift, 
progressively smaller galaxies are contributing a larger fraction of the total X-ray luminosity. 
Because X-ray selected samples trace the evolution of the overall X-ray luminosity, regardless 
of whether the contribution from different galaxy populations changes, the redshift evolution 
of an X-ray selected sample may be less dramatic than a sample which traces the evolution 
in the X-ray luminosity of a particular galaxy population. 

By selecting our red galaxies to have the same evolution-corrected, absolute A-band 
magnitudes, we are tracing the evolution in the X-ray properties of the same galaxy popu¬ 
lation as a function of redshift. We are therefore measuring an intrinsic drop in accretion 
activity of this galaxy population. As discussed in Section 5.2, more luminous unobscured 
AGN within similar host galaxies will have been excluded from our analysis due to the mod¬ 
ification in the optical color of the host galaxy due to the AGN. The redshift evolution of 
the underlying population is therefore likely to be even stronger than is observed. Since the 
overwhelming majority of luminous AGN at high redshift are hosted by luminous early-type 
galaxies (e.g., McLure et al. 1999; Kauffmann et al. 2003), we are likely to be be tracing the 
same population at lower reclshifts. The decrease in X-ray luminosity with time in our sam¬ 
ple therefore suggests we are witnessing the tailing off of the accretion activity onto SMBHs 
in earlv-type galaxies since the more powerful quasar epoch. 


6. Conclusions 

Our primary conclusions are as follows: 

• We have detected significant X-ray emission from the red galaxy population in the NDWFS 
in each of the 0.3 < z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9 redshift bins. 

• The large X-ray luminosity and hard spectral slope strongly suggest that the dominant 
X-ray emission mechanism is accretion onto a SMBH. The spectral energy distribution can 
be reproduced by invoking either a moderately obscured AGN (N H ~ 2 x 10 22 cm -2 , T=1.7) 
or a harder ADAF-like (T=0.7) spectrum. 

• The mean hard X-ray luminosity of the red galaxy population increases as a function of 
redshift. Monte Carlo simulations show that our method will correctly recover any intrinsic- 
evolution. 

• The hard X-ray emission produced by this population contributes at least 5% of the hard 
X-ray background. 

• We interpret the redshift evolution of the X-ray luminosity as a global decline in the mean 
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AGN activity of normal early-tvpe galaxies. 


One intriguing aspect of our observations is the relative hardness of the mean X-ray 
spectrum in the stacked data. We have interpreted this as resulting from either an absorbed 
AGN spectrum or a radiatively inefficient accretion flow (such as an ADAF). In the latter 
case, we envision a scenario where the accretion onto SMBHs turns progressively more ra¬ 
diatively inefficient, resulting in a harder mean spectrum and lower X-ray luminosity with 
decreasing redshift. Our results suggest that nuclear accretion is still occurring in most 
present-day early-type galaxies, but with lower efficiency and/or in a more heavily obscured 
regime. 

To distinguish between these possibilities and estimate the true accretion rate, we intend 
to extend this study to encompass the full 9 deg 2 survey area. The larger sample will contain 
an order of magnitude more galaxies. It will therefore enable us to determine the best-fit 
model to the X-ray spectral energy distribution and hence quantify the typical accretion rate 
onto SMBHs in bulge dominated galaxies and its evolution with redshift. 
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